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AB *■ Exploring the possible links between the mathematical field of 

fractional calculus and the electromagnetic theory has been one of the 
topics of our research interests. We have studied the possibility of 
bringing the tools of fractional calculus and electromagnetic theory 
together, and have explored and developed the topic of fractional paradigm 
in electromagnetic theory. Fractional calculus is a branch of mathematics 
that addresses the mathematical properties of operation of fractional 
differentiation and fractional integration operators involving derivatives 
and integrals to arbitrary non-integer orders. We have applied the tools 
of fractional calculus in various problems in electromagnetic fields and 
waves, and have obtained interesting results that highlight certain 
notable features and promising potential applications of these operators 
in electromagnetic theory. Moreover, since fractional 

derivatives/integrals are effectively the result of f ractionalization of 
differentiation and integration operators, we have investigated the notion 
of f ractionalization of some other linear operators in 

electromagnetic theory. Searching for such operator f ractionalization has 

led us to interesting solutions in radiation and scattering problems. 
CT APPROXIMATION THEORY; DIFFERENTIATION; ELECTROMAGNETIC WAVE SCATTERING; 

FRACTALS; INTEGRATION; MATHEMATICAL OPERATORS; PHYSICAL OPTICS 
ST f ractionalization methods; EM radiation problems; EM scattering problems; 

fractional calculus; electromagnetic theory; fractional paradigm; 

fractional differentiation operators; fractional integration operators; 

electromagnetic fields; electromagnetic waves; linear operators; 

physical optics approximation 
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AB Exploring the possible links between the mathematical field of 

fractional calculus and the electromagnetic theory has been one of the 
topics of our research interests. We have studied the possibility of 
bringing the tools of fractional calculus and electromagnetic theory 
together, and have explored and developed the topic of fractional paradigm 
in electromagnetic theory. Fractional calculus is a branch of mathematics 
that addresses the mathematical properties of operation of fractional 
differentiation and fractional integration operators involving derivatives 
and integrals to arbitrary non-integer orders. We have applied the tools 
of fractional calculus in various problems in electromagnetic fields and 
waves, and have obtained interesting results that highlight certain 
notable features and promising potential applications of these operators 
in electromagnetic theory. Moreover, since fractional 

derivatives/integrals are effectively the result of fractionalization of 
differentiation and integration operators, we have investigated the notion 
of fractionalization of some other linear operators in 

electromagnetic theory. Searching for such operator fractionalization has 
led us to interesting solutions in radiation and scattering problems. 



CC A4110H Electromagnetic waves: theory; A0260 Numerical approximation and 

analysis; B5210 Electromagnetic wave propagation; B0290M Numerical 

integration and differentiation; B0290F Interpolation and function 

approximation (numerical analysis) 
CT APPROXIMATION THEORY; DIFFERENTIATION; ELECTROMAGNETIC WAVE SCATTERING; 

FRACTALS; INTEGRATION; MATHEMATICAL OPERATORS; PHYSICAL OPTICS 
ST f ractionalization methods; EM radiation problems; EM scattering problems; 

fractional calculus; electromagnetic theory; fractional paradigm; 

fractional differentiation operators; fractional integration operators; 

electromagnetic fields; electromagnetic waves; linear operators; 

physical optics approximation 

L4 ANSWER 2 OF 5 INSPEC (C) 2 004 IEE on STN 

AB The family of Gowdy universes with the spatial topology of a three-torus 
is studied both classically and quantum mechanically. Starting with the 
Ashtekar formulation of Lorentzian general relativity, we introduce a 
gauge fixing procedure to remove almost all of the nonphysical degrees of 
freedom. In this way, we arrive at a reduced model that is subject only to 
one homogeneous constraint. The phase space of this model is described by 
means of a canonical set of elementary variables. These are two real, 
homogeneous variables and the Fourier coefficients for four real fields 
that are periodic in the angular coordinate which does not correspond to a 
Killing field of the Gowdy spacetimes . We also obtain the explicit 
expressions for the line element and reduced Hamiltonian. We 
then proceed to quantize the system by representing the elementary 
variables as linear operators acting on a vector space 

of analytic functionals . The inner product on that space is selected by 
imposing Lorentzian reality conditions. We find the quantum states 
annihilated by the operator that represents the homogeneous constraint of 
the model and construct with them the Hilbert space of physical 
states. Finally, we derive the general form of the quantum observables of 
the model . 

ST Gowdy universe; three- torus; spatial topology; Lorentzian general 

relativity; gauge fixing; phase space; Fourier coefficients; angular 
coordinate; Gowdy spacetimes; line element; reduced Hamiltonian; 
canonical quantization; quantum states; Hilbert space; physical 
states; quantum observables 
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AB The family of Gowdy universes with the spatial topology of a three -torus 
is studied both classically and quantum mechanically. Starting with the 
Ashtekar formulation of Lorentzian general relativity, we introduce a 
gauge fixing procedure to remove almost all of the nonphysical degrees of 
freedom. In this way, we arrive at a reduced model that is subject only to 
one homogeneous constraint. The phase space of this model is described by 
means of a canonical set of elementary variables. These are two real, 
homogeneous variables and the Fourier coefficients for four real fields 
that are periodic in the angular coordinate which does not correspond to a 
Killing field of the Gowdy spacetimes. We also obtain the explicit 
expressions for the line element and reduced Hamiltonian. We 
then proceed to quantize the system by representing the elementary 
variables as linear operators acting on a vector space 



of analytic functionals. The inner product on that space is selected by 
imposing Lorentzian reality conditions. We find the quantum states 
annihilated by the operator that represents the homogeneous constraint of 
the model and construct with them the Hilbert space of physical 
states. Finally, we derive the general form of the quantum observables of 
the model . 
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AB This paper uses a wavelet transform approach to actively image an object 
continuously distributed in range and velocity. It is shown that by 
transmitting a high resolution, i.e. large time -bandwidth product, signal 
into the environment and operating on the echo with a wavelet transform, 
an estimate of the delay-time- scale representation or wideband spreading 
function of the object can be obtained. The wideband spreading function 
(WBSF) is a characterization of the time-varying propagation/scattering 
associated with the channel being imaged. It is shown that the 
linear operator that acts on the wideband spreading 

function to form the echo is in the form of an inverse wavelet transform 
and the adjoint operator is in the form of a forward wavelet transform. 
Thus, the wavelet transform is a natural transform for the investigation 
of wideband spreading functions. By combining information extracted from 
wavelet estimates of the WBSF associated with independently located 
sensors, it is possible to estimate vectors which describe the 
physical characteristics of the object in the channel. 

Specifically, the support curve of the WBSF in the wavelet domain can be 
directly related to the projections of these vectors along the 
line of sight of each of the sensors. Therefore, it is necessary 
to obtain a wavelet estimate of the WBSF which highly resolves the support 
curve. The wavelet transform estimate is shown to be limited by the 
resolution capabilities of the auto-wavelet transform of the transmitted 
signal therefore establishing the need for high resolution transmit 
signals . 
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AB This paper uses a wavelet transform approach to actively image an object 
continuously distributed in range and velocity. It is shown that by 
transmitting a high resolution, i.e. large time-bandwidth product, signal 
into the environment and operating on the echo with a wavelet transform, 
an estimate of the delay-time-scale representation or wideband spreading 
function of the object can be obtained. The wideband spreading function 
(WBSF) is a characterization of the time-varying propagation/scattering 
associated with the channel being imaged. It is shown that the 
linear operator that acts on the wideband spreading 

function to form the echo is in the form of an inverse wavelet transform 
and the adjoint operator is in the form of a forward wavelet transform. 
Thus, the wavelet transform is a natural transform for the investigation 
of wideband spreading functions. By combining information extracted from 
wavelet estimates of the WBSF associated with independently located 
sensors, it is possible to estimate vectors which describe the 
physical characteristics of the object in the channel. 

Specifically, the support curve of the WBSF in the wavelet domain can be 
directly related to the projections of these vectors along the 
line of sight of each of the sensors. Therefore, it is necessary 
to obtain a wavelet estimate of the WBSF which highly resolves the support 
curve. The wavelet transform estimate is shown to be limited by the 
resolution capabilities of the auto-wavelet transform of the transmitted 
signal therefore establishing the need for high resolution transmit 
signals . 
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AB The authors describe a multivariable deconvolution method which 

approximates the unknown kernel by a finite series of Laguerre functions 
(which define a complete basis in L/sup 2/(0, infinity )) . The possibility 
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Eigen-values and Eigen-vectors 

Definition 4.16 A subspace JVC C V is said to be an invariant subspace of an linear operator X if 
V/ G UXf G X. 

Definition 4.17 Let T be linear operator. If f is a non-zero vector satisfying 

TJ = \f 

or some scalar A , we say that f is an eigenvector of operator T and A is the corresponding eigen- 
value. 

Note that / = 0 will always satisfy the equation Tf = A/ for an arbitrary A . Therefore, null vector 
is, by definition, excluded from being an eigen-vector. 

It is possible that for a given A there are more than one eigen-vectors satisfying the eigen- value 
equation Tf = A . Therefore, we define 

Definition 4.18 Let A be an eigen-value of an operator T. Let ^(A) denote the number of linearly 
independent eigen-vectors Tx = Ax. If =1 we say that the eigen-value A is non- 
degenerate. When l/(A) > 1, we say that the eigen-value A is degenerate and the degeneracy of the 

eigen-value A is defined to be equal to the number of linearly independent eigen-vectors with eigen- 
value A. 
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In linear algebra , the eigenvectors (from the German eigen meaning "inherent, characteristic") of a linear 
operator are non-zero vectors which, when operated on by the operator, result in a scalar multiple of themselves. 
The scalar is then called the eigenvalue associated with the eigenvector. 

In applied mathematics and physics the eigenvectors of a matrix or a differential operator often have important 
physical significance. In classical mechanics the eigenvectors of the governing equations typically correspond to 
natural modes of vibration in a body, and the eigenvalues to their frequencies. In quantum mechanics , operators 
correspond to observable variables, eigenvectors are also called eigenstates, and the eigenvalues of an 
operator represent those values of the corresponding variable that have non-zero probability of occurring. 

Examples 

Intuitively, for linear transformations of two-dimensional space R 2 , eigenvectors are thus: 

• rotation: no eigenvectors 

• reflection: eigenvectors are perpendicular and parallel to the line of symmetry, the eigenvalues are -1 and 
1 , respectively 

• scaling: all vectors are eigenvectors, and the eigenvalue is the scale factor 

• projection onto a line: eigenvectors with eigenvalue 1 are parallel to the line, eigenvectors with eigenvalue 
0 are parallel to the direction of projection 

Definition 

Formally, we define eigenvectors and eigenvalues as follows: If A : V -> V is a linear operator on some vector 
space V, v is a non-zero vector in Vend c is a scalar (possibly zero) such that 

\\mathbf{A} \\mathbf{v} = c \\mathbf{v}, 

then we say that v is an eigenvector of the operator A, and its associated eigenvalue is c. Note that if v is an 
eigenvector with eigenvalue c, then any non-zero multiple of v is also an eigenvector with eigenvalue c. In fact, 
all the eigenvectors with associated eigenvalue c, together with 0, form a subspace of V, the eigenspace for the 
eigenvalue c. 

Finding eigenvectors 

For example, consider the matrix 

A = \\begin{bmatrix} 0 & 1 & -1 \\\1 & 1 & 0 \tt-1 & 0 & 1 \\end{bmatrix} 
which represents a linear operator R 3 -> R 3 . One can check that 

A \\begin{bmatrix}1 WW 1 WW -1 Wehd{bmatrix} = 2 Wbegin{bmatrix}1 WW 1 WW -1 Wend{bmatrix} 
and therefore 2 is an eigenvalue of A and we have found a corresponding eigenvector. 

The characteristic polynomial 

An important tool for describing eigenvalues of square matrices is the characteristic poly nomial: saying that c is 
an eigenvalue of A is equivalent to stating that the system of linear equations (A - cl) x = 0 (where I is the identity 
matrix) has a non-zero solution x (namely an eigenvector), and so it is equivalent to the determinant det(A - c I) 
being zero. The function p(c) = det(A - cl) is a polynomial in c since determinants are defined as sums of 
products. This is the characteristic polynomial of A; its zeros are precisely the eigenvalues of A. If A is an n-by-n 
matrix, then its characteristic polynomial has degree n and A can therefore have at most n eigenvalues. 
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Returning to the example above, if we wanted to compute all of A's eigenvalues, we could determine the 
characteristic polynomial first: 

p(x) = \\det( A - xl) = Wdet 

Wbegin {bmatrix} 

-x & 1 & -1 \\\1 & 1-x & 0 \tt-1 & 0 & 1-x \\end{bmatrix} 

: = -x A 3 + 2x A 2 + x - 2\\ and because p(x) = -(x - 2) (x - 1) (x + 1) we see that the eigenvalues of A are 2, 1 and 
-1. 

(In practice, eigenvalues of large matrices are not computed using the characteristic polynomial. Faster and 
more numerically stable methods are available, for instance the QR decomposition .) 

Complex eigenvectors 

Note that if A is a real matrix, the characteristic polynomial will have real coefficients, but not all its roots will 
necessarily be real. The complex eigenvalues will all be associated to complex eigenvectors. 

In general, if .... v m are eigenvectors to different eigenvalues A 1( A m , then the vectors v 1( .... v m are 
necessarily linearly independent . 

The spectral theorem for symmetric matrices states that, if A is a real symmetric n-by-n matrix, then all its 
eigenvalues are real, and there exist n linearly independent eigenvectors for A which all have length 1 and are 
mutually orthogonal . 

Our example matrix from above is symmetric, and three mutually orthogonal eigenvectors of A are 

v_1 = \\begin{pmatrix} 1\\\\ 1\\\\ -1\\end{pmatrix} v_2 = \\begin{pmatrix} OWW 1\\end{pmatrix} v_3 = 
\\begin{pmatrix} 2\\\\ -1WW 1\\end{pmatrix} 

These three vectors form a basis of R 3 . With respect to this basis, the linear map represented by A takes a 
particularly simple form: every vector x in R 3 can be written uniquely as \\mathbf{x} = x_1 \\mathbf{v}_1 + x_2 
\\mathbf{v}_2 + x_3 \\mathbf{v}_3 and then we have \\mathbf{A x} = 2x_1 \\mathbf{v} 1 + x 2 \\mathbf{v} 2 - 
x_3 \\mathbf{v}_3. 

Infinite-dimensional spaces 

The concept of eigenvectors can be extended to linear operators acting on infinite-dimensional Hilbert spaces or 
Banach spaces . 

There are operators on Banach spaces which have no eigenvectors at all. For example, take the bilateral shift on 
the Hilbert space \\ell A 2(\\mathbb{Z}); it is easy to see that any potential eigenvector can't be square-summable, 
so none exist. However, any bounded linear operator on a Banach space Vdoes have non-empty spectrum. 
The spectrum Wsigma(T) of the operator T : V - V is defined as 

\\sigma(T)={\\lambda\\in\\mathbb{C}: (Wlambda 1 -T) is not invertibleW}. 

Then Wsigma(T) is a compact set of complex numbers, and it is non-empty. When T is a compact operator (and 
in particular when T is an operator between finite-dimensional spaces as above), the spectrum of T is the same 
as the set of its eigenvalues. 

The spectrum of an operator is an important property in functional analysis . 
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In mathematics, Banach spaces, named after Stefan Banach who studied them, 
are one of the central objects of study in functional analysis. Banach spaces are 
typically infinite-dimensional spaces containing functions. 

Definition 

Banach spaces are defined as complete normed vector spaces. This means that a 
Banach space is a vector space \/over the real or complex numbers with a norm 
||.|| such that every Cauchy sequence (with respect to the metric d(x, y) = ||x - y||) in 
Vhas a limit in V. 

Click the link for more information. 
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LINEAR OPERATORS ON MATRICES: PRESERVING SPECTRUM 
AND DISPLACEMENT STRUCTURE 



KENNETH R. DRIESSEL* AND WASIN SO* 



Abstract. In this paper we characterize those linear operators on general matrices that preserve 
singular values and displacement rank. We also characterize those linear operators on Hermitian 
matrices that preserve eigenvalues and displacement inertia. 

Key words, linear operator, displacement structure, Toeplitz 
AMS(MOS) subject classification. 15A04 
1. Introduction. We introduce some notation to facilitate our discussion. 



For 1 < i < m, 1 < j < n, let E % * denote the mx n matrix with 0 everywhere except 
1 at the (ij) position. Then {E^} is a basis for C mxn . We also adopt the following 



sing(>l) := the singular values of a matrix A (including multiplicity); 

eigen(^4) := the eigenvalues of a Hermitian matrix A (including multiplicity); 

rank(A) := the rank of a matrix A) 

inertia(A) := the inertia of a Hermitian matrix A. 

For A e C mxn , rank(A) = k if and only if A has exactly k nonzero singular values. 
For A e Herm{m), inertia(A) = (p, n, z) if and only if A has p positive, n negative 
and z zero eigenvalues, m = p + n + z. 

We are interested in the spectral properties of matrices that are Toeplitz or nearly 
Toeplitz. As a consequence, we are interested in linear operators that preserve these 
properties. We know of only one previous result in this direction. It is the following 
theorem due to Chu [1992]. Let E m denote the m x m exchange matrix defined by 



We write E in place of E m when m is easily determined from the context. 
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Gl{m) 
Herm(m) 
U(m) 



the set of all m x n complex matrices; 
the set of all nonsingular mx m matrices; 
the set of all m x m Hermitian matrices; 
the set of all mx m unitary matrices. 



notation. 



Em(iJ) '= 8{i,m + 1 - j). 



where 8 denotes the Kronecker delta. For example, when m = 3, 
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Theorem 1.1. Let Q be an m x m orthogonal matrix. Then the following 
conditions are equivalent: 

(*) If A is an m x m symmetric Toeplitz matrix then so is QAQ T . 
(**) The matrix Q is one of the following: 

±/, ±E, ±J', ±I'E 

where I denotes the mxm identity matrix and V := Diag{— 1, (— l) 2 , . . . , (— l) m_1 ). 

His techniques can be used to characterize nonzero linear operators on Hermitian 
matrices which preservse both eigenvalues and Toeplitz structure. We shall report 
results along this line in a future paper. 

In this report, we study the nonzero linear operators which preserve spectra and 
displacement structure. We begin by recalling the relevant definitions. Let Z m denote 
the mxm (lower) shift matrix defined by 

Z m {iJ) :=6(»,j+l). 

For example, when m = 3, 

/ 0 0 0 
Z 3 = 1 0 0 
V 0 1 0 

We write Z in place of Z m when m is easily determined from the context. Let V be 
the linear operator defined by 

V : = C mxn — > C mxn : X — ■+ X - Z m XZl. 

For A G C mxn , the displacement rank of A is defined as 

dis-rank (A) := rank (V.A). 

In the case m = n, V preserves Hermitian matrices. For A E Herm(m) 7 the displace- 
ment inertia of A is defined as 

dis-inertia (A) := inertia (V.A). 

Kailath appears to be one of the first to emphasize the importance of the displacement 
structure of matrices. We recall a few of the major results in this area in order to 
illustrate the significance of these concepts. Note that Toeplitz matrices usually have 
displacement rank 2. Hence matrices with low displacement rank are regarded as being 
"nearly Toeplitz" . The following result shows that displacement rank is preserved 
(loosely speaking) under inversion. It is from Kailath, Kung and Morf [1979]. 

Theorem 1.2. For A e Gl(m), dis-ran^A' 1 ) = dis-rank(EAE). 

The following inequality, due to Comon [1992], shows that if A has small displace- 
ment rank then so does its pseudo-inverse A+ : 

dis-rank(A + ) < 2 dis-rank(EAE). 

Note that Hermitian Toeplitz matrices usually have displacement inertia (l,l,m — 
2). Hence Hermitian matrices with low displacement inertia are regarded as being 
"nearly Toeplitz". Similar to displacement rank, displacement inertia is preserved 
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(loosely speaking) under inversion. We learned about this theorem from Tiberiu 
Constantinescu (Institute of Mathematics of the Romanian Academy of Sciences). 

Theorem 1.3. For A G Gl(m)nHerm(m), dis-inertia{A~ l ) = dis-inertia(EAE). 
Other versions of displacement structure can be defined and theorems analogous to 
the last two can often be proved too. See Chun and Kailath [1991], Heinig and Rost 
[1984]. 

The rest of this paper is organized as follows. In section 2, we shall characterize 
those linear operators on general matrices that preserve both rank and displacement 
rank. As a consequence, we obtain the characterization of those linear operators pre- 
serving singular values and displacement rank. The aim of section 3 is to characterize 
those linear operators on Hermitian matrices that preserve inertia and displacement 
inertia. We also obtain the characterization of those linear operators preserving eigen- 
values and displacement inertia. Then we have some concluding remarks in the final 
section. 

Acknowledgment. This research was done during the 1991-1992 academic year 
which the authors spent at the Institute for Mathematics and its Applications (IMA) 
at the University of Minnesota. We thank the members of the IMA for their hospital- 
ity. We thank Jack Conn (Department of Mathematics, University of Minnesota) for 
his suggestions which streamlined some of the proofs. We also thank Leiba Rodman 
(Department of Mathematics, College of William and Mary) for his comments. 

2. Preserving rank and displacement rank. In this section, we shall char- 
acterize those nonzero linear operators on C mxn that preserve both rank and dis- 
placement rank. We shall also characterize those that preserve singular values and 
displacement rank. Recall that two matrices A,B G C mxn are equivalent if there 
exist M G Glim), N G Gl(n) such that B = MAN. Note that A and B are equivalent 
if and only if ran k (A) = rank(B). The following theorem appears in Horn, Li and 
Tsing [1991]. It characterizes the linear operators preserving equivalence. 

Theorem 2.1. Let T : C mxn — ► C mXfl be a nonzero linear operator. Then the 
following conditions are equivalent: 

(*) T.A is equivalent to T.B whenever A is equivalent to B. 

(**) There exist M G G/(m), N G Gl(n) such that either, for all X G C mxr \ T.X = 
MXN or m=n and, for all X G C mxn , T.X = MX T N. 

As consequences of this result, we obtain the characterization of linear operators 
preserving rank, and those preserving singular values. 

Theorem 2.2. Let T : C mxn — ► C mxn be a nonzero linear operator. Then the 
following conditions are equivalent: 

(*) For allX G C mxn , rank(T.X) = rank(X). 

(**) There exist M G G/(m), N G Gl(n) such that either, for all X G C mxr \ T.X = 
MXN or m=n and, for all X G C mxn ,T.X = MX T N . 

Proof. (**) => (*). Direct verification. (*) => (**). If T preserves rank then it 
also preserves equivalence. Hence T has the required forms by Theorem 2.1. □ 

THEOREM 2.3. Let T : C mxn — ► C mxn be a nonzero linear operator. Then the 
following conditions are equivalent: 

(*) For all X G C mxn , sing(T.X) = sing(X). 

(**) There exist U G U{m), V G U(n) such that either, for all X G G mXn ,T.X = 
UXV or m=n and, for all X G C mxn ,T.X = UX T V. 

Proof (**) => (*). Direct verification. (*) => (**). If T preserves singular values 
then it also preserves rank. By Theorem 2.2, there exist M G Gl(m),N G Gl(n) 
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such that either T.X = MXN orm = n and T.X = MX T N. By the singular value 
decomposition, M = U1E1II2 and N = V1E2V2 where £/< G (7(ro) } V; G t/(n), £1 = 
Diag(a u . . . , a m ), E 2 = Diag(b u . . . , 6 n ). We consider the case when T.X = MXN. 
If X = Ul&iVf then T.X = £/ 1 Ei£'JE 2 V r 2 . Since sin^X) = sing(T.X), we have 
a»6j - 1. Consequently, ai = * ■ • = a m =: a, 61 = * ■ • = 6 n =: 6 and ab = 1. This 
implies that T.X = UXV where U := U1U2 and V := V1V2. The proof is similar for 
the other case. □ 

Now we characterize those linear operators on C mxn preserving both rank and 
displacement rank. For A G C, we shall use D n (X) to denote the n x n diagonal matrix 
with diagonal entries 1, A, . . . , A" -1 ; in symbols 

A* (A) -Diag{\^...,\ n ~ l ). 

Theorem 2.4. Let T : C mxn — ► C mxn be a nonzero linear operator. Then the 
following conditions are equivalent: 

(*) For allX G C mxn ,rank(T.X) = rank(X) and dis -rank (T.X) = dis-rank(X). 

(**) There exist A ^ 0 and lower triangular Toeplitz matrices M G Gl(m),N G 
Gl(n) such that either, for all X G C mxn ,T.X = D m (X)MXN T D n (X~ l ) or m=n 
and, for all X G C mxn ,T.X = D m (X)MX T N T D^X' 1 ). 

Before we prove this theorem, we need some preliminary lemmas. The first one 
is a characterization of matrices that nearly commute with the shift matrix. 

Lemma 2.5. Let B G C nxn and A ^ 0. Then the following conditions are 
equivalent: 

(*) BZ n = XZ n B. 

(**) There exists a lower triangular Toeplitz matrix L such that B = D n (X)L. 

Proof First we observe that D n (X)Z n ~ XZ n D n (X). (*) => (**). Let L := 
D n (X~ 1 )B. Then LZ n = Z n L. By comparing entries, one deduces that L is a lower 
triangular Toeplitz matrix. (**) => (*). Since L is a lower triangular Toeplitz matrix, 
there exists a polynomial p(x) such that L = p(Z n ). Hence BZ n = D n (X)LZ n = 
D n (X)p(Z n )Z n = D n (X)Z n p(Z n ) = XZ n D n (X)p(Z n ) = XZ n D n (X)L = XZ n B. □ 

Next we collect some basic results about the Kronecker product. For A G C nxn 
and B G C mxm , recall that A <g> B : C mxn — C mxn is a linear operator which may 
be defined by 

(A®B).X := BXA T . 

We prefer this "coordinate-free" definition to the usual one. (Compare Horn and 
Johnson [1991] or Graham [1981] or Lancaster and Tismenetsky [1985] ). A funda- 
mental property (which is easy to verify using this definition) is that 

(A ® B) o (C <g> D) = (AC <g> BD) 

where o denotes the composition of two operators. Moreover it can be proved that 
eigen(A ® 5) = {aify : 1 < i < n, 1 < j < m} where eigen(A) = {a* : 1 < i < n} 
and eigen(B) = {fy : 1 < j < m}. Hence tr(A ® B) = (trvl)(*r£). The next result, 
which is taken from Marcus and Moyls [1959], is a form of uniqueness for Kronecker 
product representations. 

Lemma 2.6. Let X {j Wi G C nxn and Y u Vi G C mxm . If ££ =1 X { ® Y { = 
Ei=i W< ® K*. Xi are linearly independent then each G Span{Vi, . . . , V s ). 
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Proof. Since the Xj are linearly independent, for each i there exists P» such that 
tr(PiXj) = 6(2, 7). Then, by "contraction", we have 

□ 

Corollary 2.7. Let X { G C nxn and Y- G C mxm . If{Xi} and {Yj} are linearly 
independent sets of matrices then {X{ <8> Y ; } is a linearly independent set. 
Proof. Assume J2ij QijXi <g) Yj =0. We rewrite this equation as 

^[X i ®(^a i> Y ; )] = 0. 

Use Lemma 2.6 (with all W{ = 0) to conclude that, for all i, 

0. 

i 

Since the Y,- are linearly independent, a^- = 0 for all □ 

Corollary 2.8. Ic< X G G/(n), Y G G/(m). TAen {X ® Y,Z„X ® Y,X ® 
Z m Y, Z n X <g> Z m Y] is a linearly independent set. 

Proof. Since X G Gl(n), {X,Z n X} is a linearly independent set. Similarly 
{Y,Z m Y} is a linearly independent set. Apply Corollary 2.7 to obtain the required 
result. 0 

The following result appears in Horn and Johnson [1991] and Graham [1981]. 
Lemma 2.9. Let trans := C nxn — > C nxn : X — ► X T denote the linear operator 
of taking transpose. Then trans has the following Kronecker product representation: 

n 

trans = ^ E ij <g> E ji 

where E %3 is the n x n matrix with 0 everywhere except 1 at the (i,j) position. 
Proof For X = (xij) G C nxn , note that E ji XE^ = XijE ji . Now we have 

trans.X = X T = ^ Xij E ji = ^ E ji XE ji = ^ E ij <g> E ji .X. 
ij »,;' i,j 

□ 

We adopt the convention that E ij = 0 if i > n,j > n, i < 1, or j < 1. Then it is 
easy to verify that = and Z = E^' 1 ^. With this observation, we 

are ready to prove the next lemma. 

Lemma 2.10. If P,Q } R,S e C nxn are such that 

® In ~ Z n ® Z n )o(Q®P) = (S®R)o trans o (/„ ® I n - Z n ® Z n ) 

then at least one of {P,R,S} is singular. 

Proof. Assume that P,R y S G G/(n). By Lemma 2.9, trans = . =1 ® £7\ 
Hence 

n 

(/» ® /« - Z n ® Z n )(Q ® P) = (S ® R){ E ij ® ®I n -Z n ® Z n ). 



Since 

n n 

Y (E ij ® E ji )(I n ®I n -Z n ®Z n ) = Y Eij ® Eji - EiU ~ 1} 0 Ej(i ' 1] 

ij = l i,i = l 

n 

= y Eij ® ^ ~ 0 

n 

we have 

n 

Q®P-Z n Q®Z n P = Y SEij ® " jE?W +1 ^ 1 '- 1 )). 

Note that {SE l * } is a linearly independent set. Use Lemma 2.6 to conclude that, for all 
iJiRiE^-Ed+W- 1 )) e Span{P,Z n P}An particular, RE 11 , RE 21 , R(E 12 - E 21 ) G 
Span{P,Z n P). Hence 

3 = dim Span{RE l \ RE 2 \ R(E 12 - E 21 )} < dim Span{P,Z n P} = 2. 

This is a contraction. 0 

Lemma 2.11. IfP,ReGl(m) and Q,S £ Gl(n) satisfy 

(I n ® / m - Z n ® Z m ) o (Q ® P) = (5 <g> R) o (7 n ® / m - Z n ® Z m ) 

<Aen ^ere eris* A ^ 0 and /ower triangular Toeplitz matrices N G G/(n), M G G/(rn) 
sucA that Q = DniX-^N and P = £> m (A)M. 

Proof. Note that we can rewrite the given equation as follows: 

(1) Q®P-Z n Q®Z m P = S®R-SZ n <$>RZ m . 

By Lemma 2.6, 5, 5Z n G Span{Q, Z n Q}, i.e. there exist a, /?, 7, <5 G C such that 

S=aQ + 0Z n Q and SZ n = jQ + 6Z n Q. 

Note that 5 = (a/ + PZ n )Q has rank n; hence a ^ 0. Also note that SZ n = 
(7/ + 6Z n )Q has rank n - 1; hence 7 = 0. Furthermore 0 ^ SZ n = SZ n Q and hence 
6^0. In summary, we have 

S = aQ + pZ n Q and SZ n = 8Z n Q. 

where a ^ 0 and 6^0. Similarly, we get 

R = aP + 6Z m P and RZ m = dZ m P 

where a / 0 and d ^ 0. Substituting back into the equation (1), we deduce that, by 
Corollary 2.8, /? = 0, 6 = 0 and act = 6d ~ 1. Thus 

S = aQ, SZ„ = 6Z n Q 

= aP, = dZ m P 

6 



and hence 

X~ 1 Z n Q = QZ n and XZ m P = PZ m 

where A := j = ~. By Lemma 2.5, 

Q = D n (\- l )N and P = D m (X)M 

where M, N are lower triangular Toeplitz matrices of dimension m, n respectively. □ 
We are now ready to prove Theorem 2.4. 

Proof. (**)=>(*). It is clear that T preserves rank. It remains to show that T 
preserves displacement rank. 
Case 1. T.X = D m (X)MXN T D^X' 1 ) 
By Lemma 2.5, we have 

T.X - Z m {T.X)Zl = D^MX^DniX-^-ZmD^MX^DiX-^ 

= D m (X)MXN T D n (X~ 1 ) - A" 1 D m {X)MZ m XXZl N T D{X~ l ) 
= Z) m (A)MX7V T D n (A- 1 ) - D m (A)MZ ro XZ^iV T D(A- 1 ) 
= 7) m (A)M(X - Z m X2%)N T D(\- x ). 

Hence 

cfzs-ran^T.X) = rank(TX - Z m (TX)Z^) 
= rank(X - Z m XZ%) 
= dis-rank(X). 

Case 2. T.X = £> m (A)MX T 7V T D n (A~ 1 ) 

Using an argument like Case 1, we conclude that T preserves displacement 
rank. 

(*) => (**). We assume that T is a nonzero linear operator that preserves rank 
and displacement rank. We define f : C mxn —+ C mxn by 

f := (/„ ® 7 m - Z n ® Z£) o To (7 n ® 7 m - Z n ® Z^)" 1 . 

Hence 

(/„ ® / m - Z„ ® z£) o T = f o (/„ 0 I m - Z n ® 

Since T preserves rank, by Theorem 2.2, there exist P € G/(m), Q e G/(n) such that 
either T = Q ® P or m = n and T = (Q 0 P) o trans. On the other hand, since 
T preserves displacement rank, it follows that T preserve rank. Then, by Theorem 
2.2, there exist R G G/(m), S € G/(n) such that either T = S®7£ or m = n and 
T = (S <g> R) o trans. We have 4 cases to consider. 
Case 1. T = Q ® P and f = S <8> R 

Note that (7 n ®7 m -Z n ®Z£)o(Q®P) = (S®P)o(7 n ®7 m -Z n ®Z£).Then, 
by Lemma 2.11, there exist A ^ 0 and lower triangular Toeplitz matrices 
N € G/(n),M E G/(m) such that Q = D^X-^N and P = D m (X)M. 
Consequently, T = D n (X~*)N 0 D m (A)M. 
Case 2. m = n,T = Q®P and T = (5 0 7E) o trans 

Note that (7 n 07 m -Z n 0Z£)o(Q0P) = (S® P) o trans o (7 n ® 7 m -Z„®Z£). 
Then, by Lemma 2.10, one of {P, P, 5} is singular, a contradiction. 
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Case 3. m = n, T = (Q ® P) o /rans and f = S <g> R 

Note that (7 n ®7 m -Z n ®Z£)o(Q®P)o trans = (S®#)o(7 n ®7 m -Z n ®Z£). 
Using the fact that 

(7„ ®Im-Zn® Z^) o ^rans = trans o (7„ ® I m - Z n ® Z^), 

we deduce that 

(In <8> 7 m - Z n ® Z*) °(Q®P) = (S®R)o trans o (7 n ® 7 m - Z n ® Z£). 

Then, by Lemma 2.10, one of {P, 5} is singular, a contradiction. 
Case 4. rn = n, T = (Q ® P) o trans and f = (S ® R) o trans 

Note that (7 n ® 7 m — Z n ® Z£) o (Q ® P) o ^rans = (5 ® 72) o trans o (7 n ® 
7 m — Z n ® Z^). Using the fact that 

(7 n ® 7 m - Z n ® Z£) o trans = trans o (7„ ® 7 m - Z n ® 

we deduce that 

(In ® 7 m - Z n ® z£) o (Q ® P) = (5 ® rt) o (7„ ® 7 m - Z n ® 

Then, by Lemma 2.11, there exist A ^ 0 and lower triangular Toeplitz ma- 
trices N e Gl(n),M e Gl(m) such that Q = D n (X~ 1 )N and P = D m (X)M . 
Consequently, T = D n (X~ 1 )N ® D m (X)M. □ 

□ 

Next we give the characterization of those linear operators on C nxn preserving 
both singular values and displacement rank. 

Theorem 2.12. Let T : C mxn — ► C mxn be a nonzero linear operator. Then 
the following conditions are equivalent: 

(*) For allX <E C mxn , sing(T.X) = sing(X) and dis-rank(T.X) = dis-rank(X). 

(**) There exist |A| = \n\ = 1 such that either, for all X E C mxn ,T.X = 
/iD m (A)XD n (A- 1 ) or m=n and, for all X £ C mxn ,T.X = ^D m (A)X T I? n (A- 1 ). 

Proof (**) => (*). Since |A| = \p\ = l,;zL> m (A) and T^A' 1 ) are unitary. Hence 
T preserves singular values. By Theorem 2.4, we know T also preserves displacement 
rank. (*) => (**). Since T preserves singular values, by Theorem 2.3, there exist U G 
U(m) t V eU(n) such that either T.X = UXV or m = n and T.X = UX T V. On the 
other hand, since T preserves both rank and displacement rank, by Theorem 2.4, there 
exist A ^ 0 and lower triangular Toeplitz matrices M 6 Gl(m),N € Gl(n) such that 
either T.X = Dm(X)MXN T D n (X~ l ) or m = n and T.X = D m (X)MX T N T D^X' 1 ). 
We consider the following 4 cases. 

Case 1. T.X = D m (X)MXN T D n (X~ 1 ) and T.X = UXV. 

For all X £ C™ xn ,D m (X)MXN T D n (X~ l ) = UXV. Then there exists a e C 
such that aD m (A)M = £/ and ^iV T 7) n (A~ 1 ) = V. Therefore both aD m (X)M 
and £ N T D n (X~ l ) are diagonal and soM = u7 m and N = vl n for some u,v e 
C. Moreover |A| = \a\ = \uv\ = 1. Consequently, T.X = fiD m (X)XD n (X~ 1 ) 
where /z := uv. 

Case 2. m = n, T.X = D m (A)MX7V T Z) n (A~ 1 ) and T.X = </X T V. 

For all X 6 C mxn , D m (X)M XN T D n (X~ l ) = </X T 7. Evaluating at X = 7 n , 
we get D m (X)MN T D n (X~ l ) = UV. Let W := U*D m (X)M = VD n (X)N~ T . 
Then WX = X T P^ for all X € C mxn . In particular, W commutes with every 
diagonal matix. Hence W is a scalar, and so X = X T for all X € C mxn , a 
contradiction. 
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Case 3. m = n, T.X = D m (\)MX T N T Dn^ 1 ) and T.X = J/XV. 

Using the same argument as in Case 2, we conclude that case 3 is impossible. 
Case 4. m = n, T.X = D m (X)MX T N T A^A" 1 ) and T.X = i7X T V r . 

Using the same argument as in Case 1, we conclude that T.X = ^D m (A)X T £> n (A _1 ). 

□ 

3. Preserving inertia and displacement inertia. In this section, we shall 
characterize those nonzero linear operators on Herm(n) that preserve both inertia 
and displacement inertia. We shall also characterize those preserving eigenvalues and 
displacement inertia. Recall that two matrices A,B G Herm(n) are ^-congruent if 
there exists S G Gl(n) such that B = SAS* . Note that A and B are *-congruent 
if and only if inertia(A) = inertia(B). The following theorem appears in Horn, Li 
and Tsing [1991]. It characterizes the linear operators preserving *-congruence. For 
simplicity, we write /„ as /, Z n as Z and D n (X) as D(X) in the following. 

Theorem 3.1. Let T : Herm(n) — ► Herm(n) be a nonzero linear operator. 
Then the following conditions are equivalent: 

(*) T.A is ^-congruent to T.B whenever A is ^-congruent to B. 

(**) There exists S G Gl(n) such that either, for all X G Herm(n),T.X = 
±SXS* or, for all X £ Herm(n) : T.X = ±SX T S\ 

As consequences of this result, we obtain the characterization of linear operators 
preserving inertia and those preserving eigenvalues. 

Theorem 3.2. Let T : Herm(n) — ► Herm{n) be a nonzero linear operator. 
Then the following conditions are equivalent: 

(*) For allX G Herm(n),ineriia{T.X) = inertia(X). 

(**) There exists S G Gl(n) such that either, for all X G Herm(n),T.X = SXS* 
or, for all X G Herm(n), T.X = SX T S* . 

Proof (**) =>• (*). Direct verification. (*) => (**). If T preserves inertia then 
it also preserves *-congruence. Hence, by Theorem 3.1, there exists S G Gl(n) such 
that either T.X = ±SXS m or T.X = ±SX T S*. However the cases with minus signs 
are ruled out because T preserves inertia. □ 

Theorem 3.3. Let T : Herm(n) — ► Herm(n) be a nonzero linear operator. 
Then the following conditions are equivalent: 

(*) For allX G Herm{n),eigen(T.X) = eigen(X). 

(**) There exists U G U(n) such that either, for all X G Herm(n),T.X = UXU* 
or, for all X G Herm(n), T.X = UX T U\ 

Proof. (**) => (*). Direct verification. (*) (**). If T preserves eigenvalues then 
it also preserves inertia. Hence, by Theorem 3.2, there exists S G G/(n) such that 
either T.X = SXS* or T.X = SX T S*. Since T preserves eigenvalues, eigen(I) = 
eigen(T.I) = eigen(SS*). Therefore SS* = I and so 5 G U(n). 0 

Now we characterize those linear operators on Herm(n) preserving both inertia 
and displacement inertia. 

THEOREM 3.4. Let T : Herm(n) — ► Herm(n) be a nonzero linear operator. 
Then the following conditions are equivalent: 

(*) For all X G Herm(n),inertia{T.X) = inertia(X) and dis-inertia(T.X) = 
dis-inertia(X). 

(**) There exists |A| = 1 and a lower triangular Toeplitz N G Gl{n) such that ei- 
ther, for allX G Herm(n),T.X = D(X)NXN*D(X)* or, for all X G Herm(n),T.X = 
D(X)NX T N*D(Xy. 

Proof (**) => (*). It is clear that T preserves inertia. It remains to show that T 
preserves displacement inertia. 
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Case 1. T.X = D(X)NXN*D(Xy 
By Lemma 2.5, we have 



T.X - Z(TX)Z T = D(\)NXN*D(\y - ZD{X)NXN* D(X)* Z T 

= D(X)NXN*D(Xy - X- 1 D(X)NZX(X"')- 1 Z T N*D{Xy 
= D(X)NXN*D(Xy - D(X)NZXZ T N*D(Xy 
= D(X)N(X - ZXZ T )N*D(Xy 

Hence 

dis-inertia(T.X) = inertia(T.X - Z(T.X)Z T ) 
= inertia(X - ZXZ T ) 
= dis-inertia(X) . 

Case 2. T.X = D(X)NX T N* D(Xy 

Using an argument like Case 1, we conclude that T preserves displacement 
inertia. 

(*) => (**). We assume that T is a nonzero linear operator that preserves inertia 
and displacement inertia. We define f : Herm(n) — ► Herm(n) by 

f:= (I ® I - Z ® Z T ) oT o (I ® I ~ Z ® Z 7 )- 1 . 

Hence 

(7® I - Z ®Z T )oT ^fo{I® I- Z®Z T ). 

Since T preserves inertia, by Theorem 3.2, there exist S € Gi(n) such that either 
T = 5® 5 or T* = (5® 5) o trans. On the other hand, since T preserves displacement 
inertia, it follows that f preserve inertia. Then, by Theorem 3.2, there exist R G Gl(n) 
such that either t = R® R or_T = (R<& R) o trans. We have 4 cases to consider. 
Case 1. T = 5® 5 and f= R<g>R 

Note that (I®I-Z®Z T )o(S®S) = (R®R)o(I®I-Z®Z T ). Then, by Lemma 
2.11, there exist X ^ 0 and lower triangular Toeplitz matrices N,M e Gl(n) 
such that S = D{X~ l )M and S = D(X)N. Hence M = D(\X\ 2 )N, and so 
|A| = 1. Consequently, T.X = D(X)NX(D(X)Ny for all X <E Herm(n). 
Case 2. T = S ® S and T = (7? ® R) o tran5 

Note that (7 ® 7 - Z ® Z T ) o (5 ® S) = (5 ® 72) o fran 5 o (7 ® / - Z ® Z T ). 
By Lemma 2.10, we get a contradiction. 
Case 3. T = (5 ® S) o trans and f = R® R 

Note that (7 ® 7 - Z ® Z T ) o (5 ® 5) o trans = (£ ® 7?) o (7 ® 7 - Z ® Z T ). 
Using the fact that 

(7 ® 7 - Z ® Z T ) o trans = trans o (7 ® 7 - Z ® Z T ), 

we deduce that 

(7 ® 7 - Z ® Z T ) o (5 ® 5) = (R ® A) o trans o (7 ® 7 - Z 0 Z T ), 
which leads to a contradiction by Lemma 2.10. 
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Case 4. T = (S 0 S) o trans and f = (R®R)o trans 

Note that - Z<S>Z T )o(S®S)otrans = (R®R)otranso(I®I-Z®Z T ). 
Using the fact that 

(/ 0 / - Z 0 Z T ) o ^rans = trans o (/ 0 I - Z 0 Z T ), 

we deduce that 

(/ 0 / - Z 0 Z T ) o (S 0 5) = (R 0 o (J 0 I - Z 0 Z T ) . 

Then we obtain the required result as in Case 1. 

□ 

Next we give the characterization of those linear operators on Herm{n) preserving 
both eigenvalues and displacement inertia. 

Theorem 3.5. Let T : Herm(n) — ► Herm(n) be a nonzero linear operator. 
Then the following conditions are equivalent: 

(*) For all X £ Herm(n),eigen(T.X) = eigen(X) and dis-inertia(T.X) = 
dis-inertia(X). 

(**) There exists \X\ = 1 suck that either, for all X £ Herm(n),T.X = D(X)XD(X) m 
or, for allXe Herm(n),T.X = D{X)X T £>(A)* . 

Proof (**) (*). Since |A| = 1, D(X) is unitary and hence T preserves eigen- 
values. From Theorem 3.4, it is clear that T also preserves displacement inertia. 
(*) => (**). Since T preserves eigenvalues, from Theorem 3.3, there exists U € U(n) 
such that either T.X = UXU* or T.X = UX T U\ On the other hand, since T pre- 
serves both inertia and displacement inertia, from Theorem 3.4, there exists |A| = 1 
and a lower triangular Toeplitz N £ Gl(n) such that either T.X — D(X)NXN*D(X)* 
or T.X = D(X)NX T N«D(Xy. We consider the following 4 cases. 
Case 1. T.X = D(X)NX N* D(X) m and T.X = UXU*. 

Evaluating at X = I, we get D(X)NN*D(X)* = UU* = /, i.e. D(X)N is 

unitary. Since D(X)N is lower triangular, it must be diagonal and so is N. 

This implies that N = fil because N is Toeplitz. Moreover \fi\ = 1. Finally 

T.X = D(X)XD(X)\ 
Case 2. T.X = D(X)NXN*D(Xy and T.X = UX T U\ 

For all X £ Herm(n), D(\)NXN* D(\)* = UX T U*. Evaluating at X = /, 

we find that D(X)N is unitary. Let V := U*D(X)N. Then VX = X T V for 

all X e Herm(n). This implies that V must be a scalar, and so X = X T for 

all X £ Herm(n) , a contradiction. 
Case 3. T.X = D(X)NX T N* D(X)* and T.X = UXU*. 

Using the same argument as in Case 2, we conclude that Case 3 is impossible. 
Case 4. T.X = D(X)N X T N* D(X)* and T.X = UX T U\ 

Using the same argument as in Case 1, we conclude that T.X = D(X)X T D(X)* . 

□ 

4. Concluding Remarks. There are many papers on rank preserving linear 
operators and inertia preserving linear operators: Adams [1962], Adams, Lax and 
Phillips [1965], Atkinson [1983], Atkinson and Lloyd [1980,1981], Baruch and Loewy 
[1991], Beasley [1970,1981,1983,1988], Botta [1987], Eisenbud and Harris [1988], Flan- 
ders [1962], Helton and Rodman [1985], Johnson and Pierce [1985,1986], Laffey and 
Loewy [1990], Lim [1979], Loewy [1989,1990,1991], Loewy and Pierce [1991], Marcus 
and Moyls [1959], Meshulam [1985,1989], M oore [1966], Pierce and Rodman [1988], 
Schneider [1965], Sylvester [1986], Westwick [1967,1987]. (We received this list from 
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R. Loewy. We have not looked at all these papers.) Some of them characterize lin- 
ear preservers of one particular rank class or one particular inertia class (rather than 
characterizing preservers of all rank or inertia classes as was done in Theorem 2.2 
and 3.2). These results probably make it possible to characterize linear preservers 
of one particular rank- and- displacement-rank class and linear preservers of one par- 
ticular inertia- and displacement-inertia class. Many of the results in these references 
treat rank preservers or inertia preservers over the field of real numbers (rather than 
the field of complex numbers that we used in this report). Some of the references 
even deal with more general fields of numbers. These preserver results over other 
fields probably make it possible to extend the results of this report to other fields of 
numbers. 

We mentioned earlier that there are definitions of displacement structure that 
are different than the ones we use in this report. (See Chu and Kailath [1991] and 
Heinig and Rost [1984].) There are linear preserver questions analogous to the ones 
we studied here for the other definitions. We expect that the techniques that we have 
used here can be used to easily settle such analogous questions. 

In the introduction of this report we noted that we are interested in the spectral 
properties of matrices that are Toeplitz or nearly Toeplitz. In particular, we are 
interested in sets having the following forms 

eigen~ 1 (X) 0 Toep(m) 



or 



eigen 1 (A) n dis-inertia 1 (p, n, z) 

where 



A 

eigen' 1 (X) 
Toep(m) 
dis-inertia -1 (p, n, z) 



= (Ai, . . A m ) £ R m , 

— {Ae Herm{m) : eigen(A) = A}, 

= {AeC mxm : A is Toeplitz}, 

= {A e Herrn(m) : dis-inertia(A) = (p, n, z)}. 



Now, by the spectral theorem, we have that 

eigen- 1 (X) = {QDiag(X)Q* : Q £ U(m)}. 
^From this we see that linear spectra preserving operators 

Herm(m) — ► Herm(m) : X — ► QXQ* 

for Q e U(m) can be used to move around this isospectral surface eigen~ l (X). In 
more technical language, we see that eigen' 1 (X) is the orbit of Diag(X) under the 
group action defined by 

U(m) x Herm(m) — ► Herm(m) : (Q 9 X) — ► QXQ\ 

We originally hoped that we could move around somewhat freely on the sets of the 
form eigen' 1 (X) ndis- inert ia~ 1 (p, n, z) by means of the linear preservers of such sets. 
This hope motivated our study of linear preservers. Unfortunately, our hope was too 
optimistic. Our results show that there are not enough such linear preservers. 
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09/ 676,801 

CROSS-REFERENCE TO RELATED APPLICATION 

This application is a non-provisional application of provisional application Serial No. 60/158,156 filed 
October 8, 1999. 

What is claimed is: 

1 . A method for determining configuration parameters describing a physical system, the method 
comprising the steps of 

measuring an output signal from the system in response to an input signal, 
the output signal being related to the configuration parameters 

by a linear operato r, 
and directly reconstructing each of the configuration parameters 

by applying a prescribed mathematical algorithm 

to the output signal. 

2. The method as recited in claim 1 wherein 

step of directly reconstructing includes the step of 
computing a configuration parameter function. 

3. The method as recited in claim 1 wherein 

step of directly reconstructing includes the step of 
computing a configuration kernel. 

4. The method as recited in claim 1 wherein 

step of directly reconstructing includes the step of 

computing a configuration parameter response function for 
each of the configuration parameters. 



5. A method for estimating a loop composition in terms of loop parameters representative of the loop 

composition comprising the steps of 

energizing the loop from a measurement end with an energy source, 
measuring a response signal from the loop at the measurement end, 
wherein each of the loop parameters is related to response signal by a linear operator, 



determined with reference to the linear operator, 
on the response signal. 

6. The method as recited in claim 5 wherein 

step of directly reconstructing includes 

the step of computing a loop parameter function. 

7. The method as recited in claim 5 wherein 

step of directly reconstructing includes 

the step of computing a loop kernel. 

8. The method as recited in claim 5 wherein 

step of directly reconstructing includes 
1 the step of computing a parameter response function for each of the loop parameters. 



estimating a loop composition of a subscriber loop in terms of loop parameters Xi, X2,..., Xj,..., 
X N , the loop having a frequency-domain response H(co, X 1( X 2 , . . . , Xj, , X N ) for the loop parameters, 
the method comprising the steps of 
(a) determining a range for each loop parameter Xj, 



and 



directly reconstructing each of the loop parameters 



by 



executing a prescribed mathematical algorithm, 




A method for 



(b) for each loop parameter Xi , 

generating a frequency-domain loop parameter function Fxj (co) 

wherein Fxj(co) = fc ... ... frsjXj H(co,Xi,X2,...,Xi,...,XN) dXidX2...d Xj...dXN, 

(c) generating a loop kernel k(co, p) for all loop parameters 

wherein k(co, p) = ft % 2 ... *n H(© i Xi i X 2i ,X n ) H(P,Xi,X 2 , ,X N ) dXidX 2 dX N ,Xi X 2 X N 

(d) generating a parameter response function gi (P) 

for each loop parameter from the integral relation FXj(o)) = £k(co, p) gj (p)dp, 

(e) energizing the loop from a measurement end with an energy source, 

(f) measuring a response signal Hr(co) = H(co, Xi, X 2 , . , Xj, , X N ) 

for the loop at the measurement end, and 

(g) directly determining each loop parameter Xi 

from the integral relation Xi = j H R (p) gi(P)dp . 

10. The method as recited in claim 9 wherein step (e) includes 

the step of computing the inverse of k(co, p) . 

1 1 . The method as recited in claim 9 wherein step (e) includes 

the step of computing the inverse of k(co, p) using singular value decomposition. 

12. The method as recited in claim 1 1 wherein step (f) includes 

the step of filtering noise from the response signal. 

13. A system for generating the loop composition in terms of loop parameters representative of the 
loop composition comprising 

a source of waves for energizing the loop 
from 

a measurement end, 

a detector for detecting a response signal 



from the loop at the measurement 2 end, wherein each of the loop parameters is 
related to response signal by an integral operator, and 
a reconstructor for directly reconstructing each of the loop parameters 
by executing a prescribed mathematical algorithm, 
determined with reference to the integral operator, on the response signal. 

14. The system as recited in claim 13 wherein 

reconstructor includes a processor for computing a loop parameter function. 

15. The system as recited in claim 13 wherein 

reconstructor includes a processor for computing a loop kernel. 

16. The system as recited in claim. 1 3 wherein 

reconstructor includes a processor for computing a parameter response function 
for each of the loop parameters. 



